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ABSTRACT 

Principal chiral models on a d-1 dimensional simplex are introduced and stud- 
ied analytically in the large N limit. The d — 0,2,4 and oo models are explicitly 
solved. Relationship with standard lattice models and with few-matrix systems in 
the double scaling limit are discussed. 
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The importance of understanding the large N limit of matrix- valued field mod- 
els cannot be overestimated. Not only is this the basic ingredient of the 1/N ex- 
pansion in the physically relevant case of QCD, but also the existence of large TV 
criticalities for finite values of the coupling is the starting point for the approach 
to 2-dimension quantum gravity known as the double scaling limit. Moreover we 
note that solutions of few-matrix systems may have a direct application to more 
complex systems in the context of strong-coupling expansion, since they may be 
reinterpreted as generating functionals for classes of group integrals that are re- 
quired in strong coupling calculations [1]. Unfortunately our knowledge of exact 
solutions for the large N limit of unitary matrix models is still limited. After Gross 
and Witten's solution [2] of the single-link problem, exact results were obtained 
only for the external field problem [3,4] and a few toy models (L = 3,4 chiral 
chains) [5,6]. 

In this letter we introduce a new class of lattice chiral models, whose large 
N behavior can be analyzed by solving an integral equation for the eigenvalue 
distribution of a single hermitian semi definite positive matrix. The models we are 
going to study are principal chiral models, with a global U(N) x U (N) symmetry, 
defined on a (d-l)-dimensional simplex formed by connecting in a fully symmetric 
way d vertices by ( rf 2 ) links. The partition function for such a system is 
defined to be 

/d d 
II dU i eX P N & E Tr W) + U 3 U i\- (!) 
i=l i>j=l 

Despite its apparent simplicity, this class of models includes most of the previously 
known solvable systems. As a function of the parameter d, which specifies the coor- 
dination number of lattice sites, it interpolates between the two-dimensional Gross- 
Witten model (with third order phase transition) and the infinite-dimensional mean 
field solution (showing a first order phase transition) of standard infinite volume 
lattice models. It also includes a case which, in the double scaling limit, corre- 
sponds to that of a c = 1 conformal field theory. 
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It is possible to eliminate the direct interactions among the unitary matrices 
U{ by introducing an identity in the form 



1 = 



/ dA exp -N/3 Tr [A - gti - u) 

J dAexp -N(3 Tr A4t 



(2) 



where A is a N x N complex matrix. As a consequence we obtain Z d = Z^/Zq 
where 

Z d = J ]JdUidA exp l-N(3TrAA^ + NpTrAj2 u\ + N (3Tr A^ - N 2 (3d 

i=\ \ i i 



(3) 



We now introduce the function 



F[BB ] ] = ^\og J dUexp yTr [BU ] + UB\ 



(4) 



F is a known function of the eigenvalues X{ of the hermitian semipositive definite 
matrix BB^ . More specifically in the large N limit, we know that [3,4] 



^.> = i5>+*>> I/2 -^5> 



(r + Xi) 1/2 + (r + Xj) 1/2 



-i-l < 5 > 



and there are two distinct phases: 

a) weak coupling r = 0, 

b) strong coupling Y.i( r + x i)' X ^ = 1- 
Up to irrelevant factors, it follows that 



Z d = J dB exp |-^ Tr BB] + N 2 dF(BB^) - N 2 /3dj , 



(6) 



where B replaces 2/3 A. 
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Morris [7] has shown that the angular integration can be performed in the case 
of complex matrices, and, again up to irrelevant factors, we may replace eq. (6) by 



00 , >. 

Z d = f dxi Y[(xi - Xj) exp I ~ J^Xi + N 2 dF( Xi ) - N 2 [3 \ . (7) 

In the large N limit it is legitimate to evaluate this integral by a saddle point 
method. The saddle-point equation resulting from eqs. (5) and (7) is 



\Jr + Xj 

2(3 



^_ 1 >^ (4 - d) y/r + Xi + dy/r + Xj 

3Ti 



(8) 



We introduce a new variable z% = y/r + x% and, in the large N limit, we assume 
that these eigenvalues lie in a single interval [a, b}. Denoting the large N eigenvalue 
density by p(z), eq. (8) becomes an integral equation for p(z), 



2(3 



-d 





■ j- dz'p(z') 



(d-2) 



z — z' z + z' 



(9) 



where the integration region is restricted by the condition, < a < b, with a and 
b determined dynamically. In particular, the normalization condition, 



p(z')dz' = 1, 



(10) 



must always be satisfied. Furthermore, one has the following constraint, 



/ 



^) d 4 < i. 



(u) 



with the equality holding exactly in the strong coupling region where a = y/r. 
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Let us begin by first discussing several simple cases where eq. (9) can be solved 
readily. For d — 0, the problem reduces to one with a pure gaussian interaction, 
and, by a more or less standard technique, one finds that 



M = 4^ - z\ (12) 



For d — 0, there is no weak coupling phase. Note also that, since the F-term in 
eq. (6) vanishes for d — 0, eq. (12) is obtained with a = 0. As a consequence, up 
to a constant, 

Z = exp N 2 In (3, (13) 

as expected. 

When d = 2 we obtain 



(14) 



and one may show that all results are consistent with a reinterpretation of the 
model as a Gross- Witten [2] one-plaquette model, with f3 c = \. 

When d — 3 the model can be mapped into the three-link chiral chain, which 
is known to possess a third order phase transition at f3 c — ^ [5]. 

The first non-trivial and new situation begins at d — 4. We have explicitly 
solved the d = 4 model, both in the weak and in the strong coupling phase. The 
eigenvalue density may be expressed in terms of elliptic integrals, and supplemen- 
tary conditions allow for the determination of a and b. 



Introducing the variable k(/3) = yjl — p-, in weak coupling we obtain, in terms 



■5 



of standard elliptic integrals K, II, and E, 



K(k)-^U(l-^,k) 



(15) 



with the condition An (3 = bE(k). In strong coupling we have 



^) = 2^T (z2 - fl2) " 1/2(6? - z2) " 1/a 



(b 2 -a 2 )K(k)-(z z -a z )U(l-^,k)\,\ 



2 „2\ 



(16) 

with the condition An(3 = b[E(k) — jpK(k)]. In both regimes, eq. (10) must also be 
satisfied. Closed form solutions for the constraints may be obtained at criticality: 
when (3 = (3 C = \, we get a = 0, b = n, and 



Pc(z) = ^log 



i + a/i- ^ 



7T^ 



(17) 



1-^/1-5 



Let's finally observe that a large-rf solution of eq. (9) may easily be found by 
assuming p(z) — > 5(z — z). The weak coupling solution is 



P>Pc=~ d i (18) 

and for strong coupling z — 0. Amusingly enough, this solution turns out to 
coincide with the large D = | mean field solution [8] of infinite volume principal 
chiral models on D-dimensional hypercubic lattices with the same coordination 
number as our corresponding models. 

We would like to add a few comments. Solving eq. (9) is certainly a well 
defined problem for any value of d, and in particular we expect to be able to find 
explicit solutions for simple cases, like d — 1 and d = 3. It is also possible to 
analyze eq. (9) numerically; details of our analytical and numerical techniques will 
be reported elsewhere; we only mention that for sufficiently large (3 > f3 c we can get 



z = (3d 



l + Wi- 



(3d 
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the eigenvalue distribution with desired accuracy, while near criticality convergence 
is slow: however within 1% accuracy we have evidence that f3 c = \ for all integer 
values of d [9]. It would be quite interesting to achieve more information, both 
qualitative and quantitative, on the ^-dependence of the phase transition. 

The thermodynamical quantity whose computation is easiest is the internal 
energy per unit link, w\, which may be obtained from 

b 

d(d - 1) W1 = 4^2 / dzp(z)(z 2 - r) - d - j (19) 

a 

One may then extract, in the vicinity of /3 C , the critical exponent for the specific 
heat, a. At present we know that when d — 2, a — —1, when d — 3, a — — |, when 
d — 4, a — 0, and for sufficiently large d the transition is first-order, i.e., a = 1. 

It is worth observing in this context that a more general model involving four 
unitary matrices and three couplings, interpolating between our d — 4 case and the 
4-link chiral chain, can be re-expressed as a model of two coupled complex matrices 
and admits many solvable limits, all characterized by a = 0, which corresponds to 
a c = 1 conformal field theory. 
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